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ABSTRACT 


The  development  of  a  laser  velocimeter  for  point  velocity  measure¬ 
ments  in  liquid  and  gas  flows  is  discussed.  The  principle  of  operation 
of  the  velocimeter  is  based  on  the  Doppler  effect,  and  the  velocity  meas¬ 
urements  are  accomplished  by  an  optical  homodyning  technique.  Velocity 
profiles  are  presented  for  laminar  flow  of  water  and  air  in  circular  flow 
channels.  It  was  found  that  a  correction  factor  was  required  when  meas¬ 
uring  velocities  of  fluids  having  indices  of  refraction  different  from  air. 
The  velocimeter  is  shown  to  be  capable  of  velocity  measurements  of  a 
vibrating  surface. 
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SECTION  I 
INTRODUCTION 


The  first  laser  velocimeter  (LV)  constructed  by  Yeh  and  Cummings 
(Ref.  1)  at  Columbia  University  was  used  to  measure  velocities  of  the 
order  of  10-4  m/sec.  Later  Foreman,  George,  and  Lewis  (Ref.  2) 
constructed  an  LV  which  was  capable  of  measuring  gas  velocities  from 
10" 2  to  10^  m/sec. 

The  original  purpose  of  developing  the  laser  velocimeter  at  the 
AEDC  was  to  study  flow  fields  of  low  speed  V/STOL  aircraft.  Since  the 
developmental  plan  was  established,  however,  interest  has  arisen  also 
in  the  high  speed  regime  encountered  in  shock  tubes  and  in  jet  plume 
studies,  and  at  present,  theoretical  effort  is  being  focused  on  these 
problems.  Simultaneously,  experiments  have  been  initiated  in  measur¬ 
ing  very  low  speed  boundary- layer  velocities  in  electrically  conducting 
fluids  of  MHD  flow  channels.  Because  of  the  wide  range  of  possible 
applications  of  this  instrument,  the  objective  of  this  program  is  to  define 
the  repeatability  and  resolution  of  the  instrument  and  to  extend  the  range 
from  a  few  tenths  of  a  centimeter  per  second  to  about  104  m/sec. 

The  use  of  the  LV  for  point  velocity  measurements  has  many  advan¬ 
tages  over  conventional  measuring  techniques.  The  flow  field  is  not  meas¬ 
urably  perturbed  by  the  LV,  since  only  a  minute  amount  of  the  low  power 
laser  beam  is  absorbed  by  the  fluid.  The  size  of  the  focal  point  of  a  focused 
laser  beam  can  be  as  small  as  a  few  microns  indicating  that  the  spatial 
resolution  of  the  LV  is  typically  of  the  order  of  a  few  microns.  Also,  the 
response  of  the  LV  is  linear  over  the  entire  velocity  range  of  interest  in 
normal  fluid  flow  studies.  The  range  of  velocities  capable  of  being  meas¬ 
ured  by  the  laser  technique  is  limited  by  the  frequency  response  of  the 
readout  instrumentation.  Velocity  measurements  from  a  small  fraction 
of  a  centimeter  per  second  to  103  m/sec  are  currently  possible.  The 
extension  of  the  upper  range  by  two  orders  of  magnitude  appears  theoreti¬ 
cally  possible.  The  LV  has  the  added  capability  of  measuring  rapidly 
varying  velocities  which  are  experienced  in  fully  developed  turbulence  at 
Reynolds  numbers  greater  than  10^  (Ref.  3). 


SECTION  II 

THEORY  OF  OPERATION 


The  principle  of  operation  of  the  laser  velocimeter  is  based  on  the 
Doppler  effect.  A  schematic  diagram  of  the  LV  system  is  shown  in  Fig.  1 
(Appendix),  Light  emitted  from  the  laser  is  focused  to  a  point  in  the  flowing 
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fluid  by  means  of  lens  Lj.  From  this  point,  a  small  portion  of  the  beam 
is  scattered  through  the  angle  d.  This  scattered  light  is  shifted  in  fre¬ 
quency  because  of  the  Doppler  effect  and  is  focused  by  lens  L2-  The 
frequency  shifted  light  is  then  incident  upon  a  front  surface  mirror  Mi, 
reflected  through  the  beam  splitter,  and  focused  on  the  photocathode 
surface.  The  unscattered  portion  of  the  laser  beam  is  focused  by  lens 
L3,  attenuated  by  a  neutral  density  filter,  transmitted  through  the  beam 
splitter,  reflected  at  the  front  surface  mirror  M2,  and  rotated  90  deg  by 
the  beam  splitter  so  that  it  becomes  incident  on  the  photocathode  surface. 
The  paths  of  both  beams  between  the  beam  splitter  and  the  photocathode 
must  be  coincidental.  At  the  photocathode  surface,  the  two  beams  are 
optically  homodyned  to  obtain  the  Doppler  frequency.  Hence,  the  fre¬ 
quency  detected  by  the  photomultiplier  is  proportional  to  the  velocity  at 
the  scattering  point. 


2.1  DERIVATION  OF  THE  DOPPLER  EQUATION 

Consider  a  scattering  center  fixed  in  a  reference  frame  S'  which  has 
a  velocity  V  relative  to  a  laboratory  fixed  frame  of  reference  S.  A  light 
source  of  frequency  uQ  and  wave  vector  kQ  in  S  will  appear  as  lJq  and  ko 
in  S'.  The  incident  light  is  scattered  by  the  scattering  center  in  S'.  The 
angular  frequency  of  the  scattered  light  is  identical  to  the  incident  light 
as  viewed  in  S  ",  i.  e. 

<U0'  =  Ois’  (1) 

To  obtain  the  Doppler  shift  frequency,  up,  the  frequency  of  the  scattered 
light  must  be  calculated  in  the  laboratory  system,  S,  where  uj)  is  defined 
as 


CJn'  =  COs  "  Wo  (2) 

Use  is  made  of  the  invariance  of  scaler  products  of  four  vectors, 
where  the  following  quantities  transform  as  four  vectors; 


-  (k0,  i 

c  ' 

L-' 

II 

P 

top  \ 
c  ’ 

ksfl 

11 

'tt  1 
w 

c 

-  (ks  > 

.  \ 

1 

=  (yV,  iyc) 

V 

=  (yV,  iyc)  =  (0,  iyc) 
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where 

i  =  (-1)V*  n  =  1,  2,  3,  4 

c  =  velocity  of  lighL 

and 

>  =  a  -  (4)2]  - 


Since  four-vector  scaler  products  are  invariant. 


kofl  t'  (i  -  ^0(1  "V 

(4) 

and 

=  ksfi  Vg 

(5) 

where  the  Einstein  summation  convention  has  been  employed, 
of  the  definitions  in  Eq.  (3)  into  Eqs.  (4)  and  (5)  yields 

Substitution 

k0^  i/i  —  —<0o y 

(6) 

and 

k.s^V^  =  -wsy 

(?) 

Hence, 

kopi  =  ks^  \V 

(8) 

In  three-vector  form,  Eq.  (8)  becomes 

yk0  -  V  -  ya0  =  yk3  •  V  -  y<ua 

(9) 

Equation  (2)  then  leads  to 

cod  =  (ks  „  k0)  .  y 

(10) 

Since 

|ks|  =  |kD|  -  2fn 

and 

'°-fr 

(11) 

Equation  (10)  becomes 


fD  =  sin  sin  (+  +  -§-)  (12) 

where  n  is  the  index  of  refraction  of  the  medium,  A0  is  the  wavelength  of 
the  emitted  radiation,  and  0  and  \jj  are  the  angles  defined  in  Fig.  2. 


2.2  DERIVATION  OF  THE  OPTICAL  HOMODYNE  EQUATION 

In  this  section,  a  simplified  derivation  of  the  optical  homodyne 
process  is  derived.  The  scattered  and  the  transmitted  beams  are 
assumed  to  be  aligned  so  that  they  leave  the  beam  splitter  at  the  same 
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point  and  are  incident  upon  the  same  point  of  the  photocatnode.  The  wave 
functions  for  the  transmitted  and  scattered  beams  are  approximated  by 
plane  waves  of  the  form 

=  A  e  iko*  ~  1  {13) 

$s  =  B  ei(klo+ kc)x  _  i,a)°  +  (14) 

By  the  superposition  theorem,  the  two  wave  functions  add  linearly. 

Hence, 


$  =  $T  +  $s  =  A  e lk°*  "  icUo1  +  A  e  l!k«  +  -  »(<ao+«Dh 

+  (B  -  A)  ei(k'>  +  kD)x  ~ 


05) 


The  photocathode  current  is  proportional  to  the  Poynting  flux,  time  aver¬ 
aged  over  the  optical  period,  which  in  turn  is  proportional  to  the  real 
part  of  Multiplication  of  Eq.  {15)  by  its  complex  conjugate  leads 

to 

F  =  Re  (<P  *  4>)  =  2AJ  +  2A(B  -  A)  +  (B  -  A)1  +  2AB  cos  (kDx  -  wDt)  (16) 

From  Eq.  (16),  it  is  seen  that  photocathode  tube  current  is  composed  of  a 
constant  or  d-c  term  and  an  a-c  term  that  varies  at  the  Doppler  angular 
frequency  uq.  The  laser  velocimeter  measurement  thus  consists  of 
measuring  the  frequency  of  a  signal  given  by  Eq.  (16)  and  converting  it  to 
a  velocity  by  Eq.  (12). 


SECTION  III 

EXPERIMENTAL  ARRANGEMENTS 


During  the  course  of  development  of  the  laser  velocimeter,  several 
experimental  arrangements  were  used.  First,  a  mirror  was  used  as  the 
moving  scattering  source.  The  mirror  was  attached  to  the  cone  of  a  loud 
speaker  which  was  driven  by  a  signal  generator  in  simple  harmonic 
motion.  The  experiment  was  designed  to  investigate  optical  alignment 
techniques  and  phototube  response.  Later,  velocities  of  fluid  media  were 
studied  under  various  flow  conditions. 


3.1  VIBRATING  MIRROR  EXPERIMENT 

The  optical  system  shown  in  Fig.  3  includes  a  scattering  source 
moving  in  the  simple  harmonic  motion  of  the  speaker  cone.  Thus,  the 
velocity  component  in  the  direction  of  the  scattered  light  beam  varies 
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sinusoidally*  The  time  dependence  of  the  velocity  is  given  by  Eq.  (17). 

V  =  Vm  COS  (togt)  cos  /9  (17) 

where  Ug  is  the  generator  frequency  driving  the  speaker,  and  0  is  the 
angle  between  the  velocity  and  the  scattered  beam.  From  Eqs.  (12)  and 
(16),  it  is  seen  that  the  phototube  output  should  be  a  signal  that  is  fre¬ 
quency  modulated  at  the  Doppler  frequency  at  a  modulation  rate  deter¬ 
mined  by  the  generator  frequency.  A  typical  oscillogram  of  the  photo¬ 
tube  output  for  a  speaker  frequency  of  50  Hz  is  shown  in  Fig.  4.  The  end 
points  of  the  harmonic  oscillation  are  represented  by  the  low  frequency 
portion  of  the  trace;  and  the  point  of  maximum  velocity,  the  equilibrium 
position  of  the  oscillator,  is  represented  by  the  high  frequency  portion 
of  the  trace.  All  frequencies  corresponding  to  velocities  from  zero  to 
Vxncos  £  are  present  in  the  homodyne  signal  as  shown  in  the  spectro¬ 
gram  of  the  signal  in  Fig.  5. 

Since  the  results  of  the  vibrating  mirror  studies  indicated  that  the 
LV  could  be  used  for  vibration  studies,  a  diffuse  aluminum  surface  was 
substituted  for  the  quarter  wave  mirror.  The  results  were  comparable 
to  those  obtained  in  the  mirror  experiment  with  the  exception  that  the 
amplitude  of  the  output  signal  was  lower  because  of  the  diffuse  scatter¬ 
ing  from  the  aluminum  surface. 

3.2  LIQUID  VELOCITY  MEASUREMENTS 

The  LV  was  used  to  measure  point  velocities  in  water  flow  which 
was  seeded  with  half-micron  polystyrene  latex  spheres.  The  specific 
gravity  of  the  seed  particles  is  unity;  therefore,  settling  of  these  par¬ 
ticles  in  the  chamber  did  not  occur,  A  schematic  diagram  of  the  liquid 
flow  system  is  shown  in  Fig.  6.  In  this  particular  series  of  experiments, 
the  flow  was  laminar  since  the  Reynolds  number  was  always  below  2000. 

The  optical  system  was  mounted  on  a  moveable  table  as  shown  in 
Fig.  7. 

The  system  was  traversed  so  that  the  focal  point  of  the  focusing 
lens  moved  across  the  flow  channel  diameter.  The  position  of  the  optical 
table,  and  hence  the  position  of  the  focal  point  relative  to  the  tube  walls, 
was  determined  to  within  5  x  10"^  cm  by  a  radial  micrometer. 

A  typical  oscillogram  of  the  homodyne  signal  for  laminar  flow  is 
shown  in  Fig.  8.  The  signal  was  amplified  and  filtered  to  remove  high 
frequency  noise.  A  typical  spectorgram  of  the  phototube  output  is  shown 
in  Fig.  9.  The  Doppler  signal  was  125  kHz  as  seen  from  the  upper  and 
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lower  sideband  frequencies  of  the  spectrum  analyzer.  Data  of  this  type 
were  obtained  at  0.025-cm  increments  along  the  tube  diameter.  A  plot 
of  the  velocity  profile  obtained  by  this  method  at  a  point  1,3  m  from  the 
entrance  area,  defined  in  Fig.  6,  is  shown  in  Fig.  10.  The  solid  line  is 
the  theoretically  predicted  parabolic  velocity  profile  (see  for  example 
Ref.  4)  of  a  fluid  flowing  in  a  smooth  pipe  in  laminar  flow  as  described 


where  R  is  the  tube  radius  and  is  the  centerline  velocity.  The  volu¬ 
metric  flow  rate,  W,  is  related  to  the  average  velocity,  <V>,  by 

W  =  <V>  A  (19) 

where  A  is  the  area  of  the  tube,  and<V>  is  obtained  by  integrating  Eq.  (18). 

<v>  ■  v-  B sr  /”  /  ['  -  (i)]  •***  (20) 


yielding 

<V>  -  ~f 

Substituting  Eqs.  (19)  and  (21)  into  (18)  gives 


(21) 


V  = 


(22) 


For  the  experimental  conditions  ip  -  90  deg,  6  ~  10.  3  deg,  n  =  1.  33 
and  A.o -6328  A,  the  centerline  Doppler  frequency  was  measured  to  be 
24.  5  kHz  which  corresponds  to  a  Vm  of  8.  68  cm/ sec.  Vm  calculated 
from  Eq.  (22)  is  8.  52  cm/ sec.  The  error  of  1.  7  percent  is  well  within 
the  accuracy  of  the  flow  rate  measurement. 


When  comparing  the  theoretical  velocity  profiles  with  experimentally 
measured  values,  a  correction  factor  for  the  distance  traversed  is 
required  because  of  the  index  of  refraction  of  the  fluid  media.  For  fluid 
media  with  indices  of  refraction  greater  than  that  of  air,  the  distance 
traversed  by  the  focal  point  of  the  laser  beam  in  the  fluid  appears  con¬ 
tracted  when  compared  to  the  actual  distance  that  the  optics  are  moved. 
This  is  seen  by  considering  the  geometry  shown  in- Fig.  11. 


The  distance  y  traversed  by  the  optical  table  is 


where 


(23) 

(24) 
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From  Snell’s  law 
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n 


Sin  g 
Sin  6 


(25) 


Thus  the  ratio  of  the  table  movement  distance  to  the  actual  traverse 
distance  in  the  medium,  obtained  by  combining  Eqs.  (23),  (24),  and 
(25),  is 


V 

D 


(26) 


For  air  flow  n  is  equal  to  unity,  and  from  Eq.  (25)  it  can  be  shown  that 

lim  =1  (27) 

n-i  \t>/ 

A  plot  of  ^  as  a  function  of  the  scattering  angle,  a,  (for  water  flow 
n  =  1.  33)  is  shown  in  Fig.  12.  This  correction  factor  was  applied  to  the 
data  in  Fig.  10. 


3.3  GAS  VELOCITY  MEASUREMENTS 

The  LV  was  used  to  measure  point  velocities  in  smoke  seeded  air 
under  laminar  flow  conditions.  The  flow  system  used  in  the  experiments 
is  shown  in  Fig.  13.  A  typical  spectrogram  of  the  homodyne  signal  is 
shown  in  Fig.  14,  and  a  plot  of  the  measured  velocity  profile  and  the 
velocity  profile  calculated  from  Eq.  (22)  is  shown  in  Fig.  15.  Again  the 
experimental  agreement  is  within  the  accuracy  of  the  measured  flow  rate. 


SECTION  IV 
CONCLUSION 


The  velocity  of  a  vibrating  mirror  and  point  velocity  measurements  in 
water  and  gas  laminar  flows  have  been  obtained  using  the  LV.  The  data 
obtained  for  the  fluid  flows  were  experimentally  verified  by  independent 
means.  The  velocity  profiles  of  the  two  experiments  were  found  to  deviate 
by  less  than  5  percent  from  the  theoretical  curves,  which  is  well  within 
the  experimental  error  involved  with  measuring  the  angles  9  and  0  and 
the  volumetric  flow  rate  W.  Since  the  quantitues  n  and  \0  Eq*  (12) 
are  accurately  known,  the  accuracy  of  the  LV  is  determined  by  the  accu¬ 
racy  to  which  the  angles  6  and  ^  can  be  measured. 

Efforts  are  underway  to  measure  velocity  distributions  at  points  in 
a  fluid  in  fully  developed  turbulence.  Coupled  with  this  effort,  a  direct 
readout  system  for  obtaining  real  time  velocity  measurements  will  also 
be  developed.  Plans  call  for  extending  the  range  of  the  LV  velocity  meas¬ 
urements  to  approximately  104  m/sec. 
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Fig.  1  Schematic  Diagram  of  a  Typical  Loser  Velocimeter  System 
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Fig.  2  Angle  Definitions  for  Doppler  Effect  Calculation 
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Fig.  3  Optical  System  far  the  Vi 
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Fig.  4  Oscillogram  of  the  Vibrating  Mirror  Signal 
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Fig.  5  Spectrogram  of  the  Vibrating  Mirror  Signal 
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Fig.  6  Schematic  Diagram  of  Liquid  Flow  System 
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Fig.  8  Typicol  Spectrogram  from  the  Water  Velocity  Measurements 
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Fig.  9  Typicol  Photomultiplier  Output  Signal 
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Fig.  10  Velocity  Profile  of  Water  in  o  Circular  Pipe 
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Fig,  11  Geometry  for  Correction  Foctor  Calculation 
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Fig.  13  Flaw  System  for  Velocity  Measurements  in  Air 
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Fig.  14  Typical  Spectrogram  from  the  Air  Velocity  Measurements 
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Fig.  15  Velocity  Profile  of  Air  in  a  Circular  Pipe 
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